Bose-Einstein condensation in a magnetic double-well potential 



T.G. Tiecke, M. Kemmann, Ch. Buggle, I. Shvarchuck, W. von Klitzing, and J.T.M. Walraven 
FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands 

(Dated: February 1, 2008) 

We present the first experimental realisation of Bose-Einstein condensation in a purely magnetic 
double-well potential. This has been realised by combining a static Ioffe-Pritchard trap with a time 
orbiting potential (TOP). The double trap can be rapidly switched to a single harmonic trap of iden- 
tical oscillation frequencies thus accelerating the two condensates towards each other. Furthermore, 
we show that time averaged potentials can be used as a means to control the radial confinement of 
the atoms. Manipulation of the radial confinement allows vortices and radial quadrupole oscillations 
to be excited. 
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I. INTRODUCTION 

Bose-Einstein condensates in dilute vapours (BEC) 
have been studied using magnetic jl], ||, [| and optical 
||, H potentials. Magnetic traps can be divided into two 
classes: static and dynamic traps. Bose-Einstein conden- 
sation has first been observed in a dynamic trap using 
a time orbiting potential (TOP) |jj and shortly there- 
after in a static Ioffe-Pritchard trap A TOP trap 
uses a three dimensional quadrupole field to confine the 
atoms. In order to avoid Majorana losses near the centre 
of the trap a rotating homogeneous field shifts the field 
minimum onto a circle around the trapped sample. At 
an appropriate rotation frequency the trapping poten- 
tial is simply the time average of the orbiting potential. 
The Ioffe-Pritchard trap on the other hand is designed to 
have a static, non-zero field minimum: in two directions 
(x,y) the atoms are confined by a quadrupole field. On 
the z-axis the radial field component is zero and the axial 
component is a parabola B z = Bq + f3z 2 with a curvature 
(3 and an offset B . 

In this paper we explore the possibilities offered by 
time averaged potentials (TAP) in Ioffe-Pritchard traps. 
We distinguish two cases Bq > and Bq < 0. In the 
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FIG. 1: Absorption image of two Bose-Einstein condensates 
created in a double TOP trap. The condensates have been 
accelerated towards each other by switching to a single trap. 
The picture is taken after 3 ms of free expansion. The figure 
on the right is an optical density profile through the centre of 
both condensates. The solid line is a fit to the experimental 
data. 



first case B z is always positive and the axial confining 
potential is a parabola. We can influence the radial and 
axial confinement of the atoms by applying a linear TAP- 
field, i.e. a homogeneous modulation field orthogonal to 
the z-axis. This allows us to introduce a radial ellipticity 
in the otherwise axially symmetric trap, which can be 
exploited to drive vortices or quadrupole oscillations. In 
the second case (B < 0), B z crosses zero at two points 
on the z-axis creating an axial double well potential. We 
can create a double TOP trap by adding a circular TAP- 
field orthogonally to the z-axis pi 0]. As an example 
we demonstrate Bose-Einstein condensation in a double 
TOP trap (see figure @). 



II. TRAP GEOMETRY 

We explored the possibilities offered by TAP fields in a 
Ioffe-Pritchard coil configuration using an efficient semi- 
analytic model of the trap. We calculate the magnetic 
field for the detailed geometry of each coil separately over 
a three-dimensional grid around the trap centre. Two 
different regions are chosen: one for the cold atom cloud 
and the Bose-Einstein condensate (10 x 10 x 100 /.im) and 
one for the hot thermal cloud (1 x 1 x 10 mm). On this 
grid a three-dimensional fourth-order polynomial is fitted 
for each of the coils. The final field is calculated as the 
sum over the polynomials with the electrical currents of 
the individual coils as parameters. The frequencies and 
nonlinearities of the trap are then simply the appropriate 
coefficients of the final polynomial. Using these polyno- 
mials any field possible with our coil configuration can be 
rapidly calculated. For selected cases we also derive sim- 
plified analytical expressions to emphasise the principles 
involved. 

In order to elucidate the TAP principle we write the 
potential U(x, y, z) of the axially symmetric (lu p = lo x = 
UJ V ) Ioffe-Pritchard trap in the approximation of a elon- 
gated trapping potential (u> p 3> w z ): 



U(x,y,z) = ^]la 2 {x 2 +y 2 ) + ( B + -f3z 2 \ (1) 
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where a is the radial gradient, p the magnitude of the 
magnetic moment of the atoms and m the atomic mass. 
For the standard Ioffe-Pritchard trap (Bo > 0) the har- 
monic frequencies are 



U! a = \ — 



p or 



rn Bq 



m 



(2) 



A. Magnetic double- well potential (Bo < 0) 

For a negative offset the modulus of the field becomes 
zero at two points on the z-axis creating two three- 
dimensional quadrupole traps. The distance Az between 
the two trap centres can be controlled via B , 



Az = 2^/2\B a \/f3 



(3) 



Inversely, we can conveniently determine Bq from a mea- 
surement of Az. 

To eliminate depolarisation near the trap minima we 
use a circular TAP, creating a double TOP configuration 
H fij. A homogeneous field of amplitude B m orthogo- 
nal to the 2-axis displaces the magnetic field zero by a 
distance 



trap to a static Ioffe-Pritchard trap. This is done by 
removing the TAP field and switching from Bo < to 
B' > 0, changing neither a nor (3. The condensates then 
accelerate towards each other and collide in the origin at 
a relative kinetic energy of -Ekm = 4/x|I?o|. In a harmonic 
trap the shape of a cloud is completely decoupled from 
its centre-of-mass motion. Therefore, if the trapping fre- 
quencies in the single and double traps are equal, the gas 
clouds move towards each other without changing shape. 
It follows from equations [2] and || that the axial and ra- 
dial frequencies of the double and single traps are the 
same if the offset and modulation fields obey 



B — 2B m — —4Bq 



(7) 



In our trap we can easily explore collision energies up to 
1.3 mK corresponding to 7000 recoil energies at 780 nm. 

Let us turn to the question of gravity. If the trap is 
slightly tilted gravity will affect the relative depth of the 
two potential minima and thus the relative number of 
atoms in the two traps. As the distance between the 
traps increases the difference in atom numbers continues 
to grow. If the traps were to remain coupled all atoms 
would eventually collect in the lower trap. Therefore the 
two traps have to separated sufficiently fast. 
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This expression is conveniently used to calibrate B m . 

Rotating the modulation field around the z-axis (with 
an angular modulation frequency u> m ) moves the zero 
along a so-called 'circle of death'. If the modulation fre- 
quency is large compared to the oscillation frequency of 
the atoms in the trap but slow compared to the Larmor 
frequency (uj PyZ -C w m <C pB m /h) the atoms will see the 
time-average of the modulated potential Q] : 



B. Radial Ellipticity (B > 0) 

For a positive offset we have a standard Ioffe-Pritchard 
potential. If we then apply a linear TAP in the IE- 
direction modulating at an angular frequency uo m we can 
write the time-averaged potential as 



U T Ap(x,y 1 z) 
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U(p rn smuj m t, y, z)dt . (8) 
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The harmonic trapping frequencies for both minima of 
the double-well potential are 
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Note that ui z = for Bo = and B m > 0. For this case 
the axial confinement is governed purely by higher order 
terms. 

The two traps are separated by a potential barrier of 
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where the first term in the brackets corresponds to the 
field at the origin and B m to the field minimum of each 
of the wells. 

One of the interesting aspects of this double trap is 
the possibility of rapidly switching from the double TOP 



The radial harmonic trap frequencies u> y ) can be 
expressed in terms of the ratio b = B m /B : 
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where K'(&) 
given by 



± Jo (1 + 6 2 sin 2 0)- 1/2 # and E'(b) = 
n 2 ^) 1 / 2 dej) B. For b <C 1 the ellipticity is 
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For b = 1 we have e = —0.16 . 

Radial ellipticity can be used for example to excite 
quadrupole oscillations or to create vortices in a Bose- 
Einstein condensate. To drive the radial quadrupole os- 
cillation the elipticity of the trap has to be modulated for 
example by driving the fields as B x (t) — B m shnv m t (1 + 
sin ilt)/2 and B y (t) — B m cosu m t (1 — sin £!i)/2, where f2 
is matched to the frequency of the quadrupole oscillation. 
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Vortices can be created by rotating an elliptic trap- 
ping potential at an angular frequency f2. This has been 
demonstrated using rotating laser beams in a mag netic 
trap j|, and in an asymmetric TOP trap jn], pj. We 
can do this by rotating the linear TAP-field as B x (t) — 
BmSmumt sin fit and B y (t) = B m $muj m t cos fit. 



III. BEC IN A DOUBLE TOP 

As an example of the use of time-averaged potentials 
in combination with a Ioffe-Pritchard coil configuration 
we demonstrate Bose-Einstein condensation in a double 
TOP trap. 

Our Ioffe-Pritchard trap is shown in figure ^ and is 
described in detail in reference juj). The trap consists 
of four Ioffe coils producing the radial quadrupolc field. 
The field of the pinch coils provides the axial parabolic 
potential. The compensation coils are designed to com- 
pensate the field of the pinch coils at the origin. Both 
sets of coils are connected in series in order to reduce 
field noise. The offset B is controlled by applying a 
dedicated homogeneous field (coils not shown in figure 
||) . In order to be able to create the modulation fields we 
have attached modulation coils to the outer faces of the 
Ioffe coils. They consist of PCB boards with one layer 
of 35 /im copper on either side, connected in series. In 
order to maximise the trap stability facing modulation 
coils were driven in series. 

The experimental path towards BEC into the mag- 
netic double-well is as follows. We optically pump about 
6 x 10 9 atoms from a 87 Rb MOT into the \5 2 S 1/2 ,F = 
2, mp — 2} state. These are loaded into a match- 
ing, horizontal Ioffe-Pritchard trap (u> p = 2tt x 8 Hz, 
uj z = 2tt x 7 Hz and B Q = 37 G). The circular TAP-field 



of B m — 0.68(3) G rotating at a frequency of 7.0 kHz 
remains constant throughout the experiment 1 14 . The 
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FIG. 2: The coil configuration of our Ioffe-Pritchard trap. 
The arrows indicate the direction of the currents in the main 
coils. Modulation coils (not shown) are attached to the out- 
side faces of the Ioffe coils. 



transfer efficiency is about 50% and the temperature in 
the trap is about 70 /.tK. We then compress the trap 
to lu p = 2tt x 390(10) Hz and uj z = 2tt x 14.8(4) Hz by 
ramping the currents in the main coils from 50 A up 
to their final value of 400 A and reducing the offset to 
Bq = 0.4 G. At this current we have a — 352(1) Gem -1 
and = 266(2) G enr 2 @ . 

Rapid forced evaporation using a radio-frequency 
sweep-down to an intermediate value of 433 kHz above 
the bottom of the trap cools the sample to a temper- 
ature of about 8fiK. At this temperature 90% of the 
atoms lie within the future circle of death. We then ramp 
the offset-field Bo linearly down at a rate of — 0.5Gs _1 
eventually splitting the cloud. 

Figure || shows a profile through a thermal cloud of 
10 6 atoms at a temperature of 0.66 fiK during the early 
stages of the splitting process at B = —53(1) mG. The 
temperature and atom number was determined by time- 
of-flight imaging. The solid curve is calculated for the 
measured temperature using the semi-analytical model 
adjusting the maximum optical density and the tilt angle 
(0.4 degrees) to fit the data. 

After the cloud is fully split we fix the offset at Bq = 
— 0.63(1) G and condense the sample by lowering the 
radio-frequency to 23 kHz above the bottom of the trap. 
Figure ^ shows a profile of an absorption image of one of 
the two Bose-Einstein condensates after a free expansion 
of 4 ms. The solid line is a fit to the Thomas- Fermi dis- 
tribution in a harmonic trap. The condensate contains 
4 x 10 4 atoms. 
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FIG. 3: An optical density profile along the z-axis of the ther- 
mal cloud during the early stages of separation after a free 
expansion of 4 ms. The cloud contains 10 6 atoms at a tem- 
perature of 0.66 /iK as determined from the radial expansion. 
The solid line is a prediction from the model with the an- 
gle (0.4 degrees) and maximal optical density as adjustable 
parameters. The inset shows the absorption image. 
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FIG. 4: An optical density profile along the z-axis of one of 
the two Bose-Einstein condensates after 4 ms of free expan- 
sion. The solid line is a fit of a Thomas-Fermi profile to the 
experimental points. The Bose-Einstein condensate contains 
4 x 10 4 atoms. The inset shows the absorption image. 




100 200 300 400 



Holding time [ms] 



FIG. 5: An axial centre-of-mass oscillation of a Bose-Einstein 
condensate in a double-well trap. The solid line corresponds 
to a fit of a sine function to the data (lo z = 2n x 26.8 Hz). 



IV. CONCLUSIONS AND OUTLOOK 



We excited a centre-of-mass oscillation of a Bose- 
Einstein condensate along the z-axis by jumping the off- 
set B from -0.43(1) G to -0.63(1) G thus shifting the 
traps by 115(1) /im outwards. The condensates started 
to oscillate around the new trap centre (dots in figure 
|[). The solid line is a fit of a sine wave to the data 
[u) z = 2-7T x 26.8 Hz). From our semi-analytical model in- 
cluding an anharmonic shift of 0.8 Hz we find a frequency 
of 27.5(8) Hz. Hence we have agreement with the data 
within the experimental error. 



We explored the TAP principle as a powerful tool to 
manipulate Ioffe-Pritchard traps. For a positive Bq a ra- 
dial ellipticity can be induced and used to excite vortices 
and quadrupole oscillations. For a negative Bq double 
well potentials are obtained. These can readily be con- 
verted to a single trap and thus used in collision exper- 
iments with ultra-cold gas clouds. As an example we 
used a circular TAP field to demonstrate Bose-Einstein 
condensation in a double TOP trap. 

This work is part of the Cold Atoms program of FOM. 
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